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We study the scattering lengths of charmonia {J/tp and rjc) with light hadrons (tt, p and A^) by 
the quenched lattice QCD simulations on 24^ x 48, 32^ x 48 and 48^ x 48 lattices with the lattice 
spacing a ~ 0.068 fm. The scattering length is extracted by using the Liischer's phase-shift formula 
together with the measurement of the energy shift AE of two hadrons on the lattice. We find 
that there exist attractive interactions in all channels, J/tp{ric)-'K, J/tp{r]c)-p and J/tp{r]c)-N: The 
s-wave J/i/j-TT {rjc-Tr) scattering length is determined as 0.0119±0.0039 fm (0.0113±0.0035 fm) and 
the corresponding elastic cross section at the threshold becomes 0.0181qqjq mb (O.OIBIqqos mb). 
Also, the J/tp-N (rjc-N) spin-averaged scattering length is 0.71±0.48 fm (0.70±0.66 fm), which is 
at least an order of magnitude larger than the charmonium-pion scattering length. The volume 
dependence of the energy shifts is also investigated to check the expected behavior of AE at 

a large spatial size L. 

PACS numbers; ll.15.Ha, 12.38.-t 12.38. Gc 



I. INTRODUCTION 



Properties of single hadrons have been studied in quantum chromodynamics (QCD) by using various techniques 
such as the QCD sum rules and lattice QCD simulations. On the other hand, the interaction between color-singlet 
hadrons is not fully explored yet because of their complex nature originating from quark exchanges and multiple gluon 
exchanges. 

In this paper, we study the low-energy elastic scattering of charmonia {J/ip and rjc) with iso-nonsinglet hadrons 
composed of up and down quarks (tt, p and N) on the basis of the quenched lattice QCD simulations. Since the 
valence quarks in each hadron in the initial state stay in the same hadron in the final state, the process we consider is 
much simpler than the interactions between light- hadrons, while it is still non-trivial in the sense that non-perturbative 
gluon exchanges play essential roles. Furthermore, such interaction has a direct relation to the physics of charmonium- 
nucleus (A) bound states 0, 0, 0, 0, Q and an indirect relation to the the elastic and inelastic charmonium-hadron 
interactions at high energies la S IS ■ 

The charmonium-iV interaction at low energies has been discussed in the framework of operator product expansion 
0,0,0,0,0,0] and in some hadronic model 8']. For example, Hayashigaki has shown that the s-wave J/tp-N 
scattering length Oq^^'^ is about 0.1 fm by using QCD sum rules while Brodsky and Miller found that it is 
about 0.25 fm from the gluonic van der Waals interaction 4J. Recently, Sibirtsev and Voloshin suggested that a lower 
bound of the Oq^^'^ is as large as 0.37 fm from the multipole expansion analysis of chromo-polarizability |13|. In the 
above-mentioned calculations, the elastic cross section at the threshold reads 1.3 mb, 7.9 mb and 17 mb, respectively. 
Brodsky, Schmidt and de Teramond introduced an effective charmonium-nucleon potential V{r) = —ae~^'^/r with 
the parameters a = 0.42 ~ 0.59 and p, = 6.0 GeV estimated by a phenomenological model of Ponieron interactions jj. 
They argued that the J/ijj-A bound system may be realized for mass number A> 3 if the attraction is sufficiently large, 
which was later confirmed by Wasson, who solved the Schrodinger equation for the charmonium- nucleus system pf. 

Unlike the case of J/ip-N, the interaction of the charmonium with the pion has a special feature due to the Nambu- 
Goldstone nature of the pion. The current algebra shows that the s-wave scattering length of a heavy hadron (H) 
with the pion is given by the well-known formula IB] : 

where /jr (Ih) is the isospin vector of tt (H) and Itv ■ Ih = \ ~^ 1) ~ Ih{Ih + 1) — 2] with I being the total 
isospin of the ■n-H system. For the tt-TV case, this is known as the Tomozawa- Weinberg relation. Eq. ^ indicates 
that the soft pion decouples from any hadrons in the chiral limit (Mtt 0). More importantly in our context, the 
leading term in Eq. iQJ vanishes even if Af^r is finite as long as H is iso-scalar (such as J/ip and rjc): Namely, the 
low-energy pion interaction with the charmonium states is very weak of 0{M^). An alternative estimate based on 
the color-dipole description of heavy quarkonia shows a small but attractive interaction between the charmonium and 



2 



the pion [lOl llllIT^ . Here, the contribution stems from the QCD trace anomaly and the effect starts from 0{M^) as 
consistent with the current algebra analysis [l7l| . 

The method we employ for extracting the charmonium-hadron scattering lengths is the quenched lattice QCD 
simulations together with the phase-shift formula by Liischer "iS] . The basic idea is to put two hadrons in a finite box 
with a spatial size L and to measure the energy of interacting hadrons relative to the energy of the non-interacting 
hadrons. Such an energy shift AE can be translated into the scattering phase shift at low energies and hence to 
the s-wave scattering length in the limit of zero relative momentum. Applications of Liischer's formula to numerical 
simulations in lattice QCD have been previously done for tt-tt, tt- N and N-N systems in Refs. ,19, 2_0„ 21. 22, 2^ 
and also for a hypothetical bound state system in Ref . . Our lattice data show that the low-energy interaction of 
the charmonia with light hadrons is always attractive. Among others, the charmonium-7r interaction is quite small 
in accordance with the current algebra result, while the charmonium-p and charmonium-A^ interactions are likely to 
be large and may even support the charmonium-nuclcus bound states according to the phenomcnological analyses 
mentioned above. 

This paper is organized as follows. In Sec. II, we recapitulate the method to calculate the s-wave scattering 
phase shift in lattice QCD and discuss its application of the Liischer's phase-shift formula to the case of attractive 
interactions. In Sec. Ill, details of our Monte Carlo simulations and numerical results of J/'(/;-hadron and r/c-hadron 
interactions at low energy are given. Physical implications of our results are also discussed. Finally, Sec. IV is devoted 
to summary and concluding remarks. In Appendix A, some details of the spin projections for J/ip-p and J/ip-N 
systems are described. 



II. TWO HADRONS IN A FINITE BOX 



It has been shown by Liischer that the s-wave scattering phase shift is related to the energy shift AE in the 
total energy of two hadrons in a finite box \lS\. To measure the total energy of two hadrons {hi and /i2) in the 
center-of-mass frame, we define the following four-point correlation function: 



(2) 



where each hadron is projected onto the zero momentum state by the summation over all spatial coordinates x, 
i.e. 0^{t) — J2x^^i^'^)- avoid the Fierz re-arrangement of two-hadron operators, we choose t4 = + 1 and 
t2 — ti + 1. Let us first consider here the scattering processes, J/ip-ir, J/tp-p and J/ip-N (Cases for 77c-7r, r]c-p and 
ric-N will be discussed in Sec. IIII Accordingly, the four-point correlation functions are defined by 



G 
G 



'(i,tsrc) = ^af(t+l)Of'^(t)(0^(isrc + l)0/^'^(W)^^, 

Here we use the conventional interpolating operators, 

Oy^^ix.t) = c„(x,t)7^Ca(x,t), 

0'^(x,t) = Ua(f,t)75da(x,t), 

0^{x,t) = eabc u^{x,t)C-f^db{x,t) Uc{x,t), 



(3) 
(4) 
(5) 



(6) 
(7) 

(8) 

(9) 



for J/i/j, TT, p and A'', respectively. C is the charge conjugation matrix, C = 7472- Also, a,b and c are color indices, 
fc and I are spatial Lorentz indices, and u, d and c are the up, down and charm quark fields. For the s-wave 
J/V'-TT scattering, the total spin is restricted to be 1. Then we simply choose diagonal correlation averaged over the 
spatial indices, i G'H^'^ {t, tsrc), to extract the scattering length. On the other hand, for s-wave J/tp-p and J/ip-N 
scatterings, different total spin states are allowed: Spin-0, 1 and 2 states for J/'ip-p and spin-1/2 and 3/2 states 
for J/ip-N . Therefore, we need appropriate spin projections to disentangle each spin contribution from four-point 
correlation functions. Details of such spin projection are described in Appendix A. 
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The hadronic two-point functions are defined as Gh{t,ts,c) = (0''(i)0'''f(isrc))- For the vector mesons such as J/ip 
and /9, we take an average over the spatial Lorentz indices as ^J2i=i 23 ^*-^i'(^)^f^(^src)^ so as to obtain possible 
reduction of statistical errors. Note that we have neglected the disconnected diagrams (such as self-annihilation of J/ifi) 
in evaluating the two-point function of J/4' and the four-point functions with J/iJj in our simulations. Contributions 
from the disconnected dia,grams in the vector channel are known to be negligibly small for strange and charm quarks 
in numerical simulations |25l l26j| in accordance with the mechanism of Okubo-Zweig-Iizuka (OZI) suppression. 

Let us now define a ratio R'^/^~^{t) to extract the total energy of two hadrons {Ej/^_i^) relative to the total energy 
of individual hadron {Mj/^, and Mh)'. 

Rj/,-n{t)- ^ r.fTrf;lTl -.n ^ ^M-^E-t), (lo) 

Gj/ip(t, tsic)Uh(t + i, Esrc + tj *>tBrc 

where 

AE = Ej/^_h - [Mj,^, + Mh). (11) 
We can introduce their relative momentum k outside the interaction range through 



where k should vanish as 1/L with increasing L if there is no bound state in this channel. 

Assuming that the interaction range R is smaller than a half of the lattice size, R < L/2, the s-wave phase shift in 
a finite box, So{k), may be written as |l8j | 

^^^Ti^ - kcotSoik) = - + 0{ei (13) 
Lit'-''' ao 

where q{= (|^)^) takes a non-integer value due to the two-particle interaction. The function 2^00(119) is an analytic 
continuation of the generalized zeta function, Zoo{s,q) = X^nez^ ("^^ ^ 9)": from the region s > 3/2 to s = 1. 
The s-wave scattering length qq is defined through the small k limit of the above formula. 

If ao/L is sufficiently small, then one can make a Taylor expansion of the phase-shift formula (|13|l around ^ 
and obtain 

^E = -'-^(l + c,^+c, r-^Y) + 0{L^) (14) 



/iL3 \ ' ^ L ' V L 

with ci = —2.837297 and C2 — 6.375183. /i denotes the reduced mass of two hadrons, /i — Mj/^AIii/{Mj/^ + M^). 
For AE > 0, Eq. (|14|) with an expansion up to 0{L~^) and that up to 0{L~^) have a real and negative solution for 
ao- On the other hand, for AE < 0, the expansion up to 0(L~'*) gives no real solution for 

although the expansion up to 0{L^^) always has a real solution. Therefore, we can use Eq. IjlSI) as a necessary 
condition to test the convergence of the large-L expansion. If this condition is not satisfied or marginally satisfied, 
we need to use the full expression of ifT^ to extract ao from AE. Eqs. (|13I15() are the basic formulas to be utilized 
in the following sections. 

III. NUMERICAL RESULTS AND DISCUSSIONS 

We have performed simulations in quenched QCD with the single plaquette gauge action 

Sg^l3Y,h-lRe[TTuA, (16) 

where /3 = 6/(7^, Up — U^i,{n) and ~ S^<iyi ^''^'^ the Wilson fermion action 

Sf = X!^'^'^) 1 ~ [t^-*- ~ lt^Wt^{n)5n+i_i,ni + (1 + lti)Ul{m)5n-i_ijn \q{.m), (17) 
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where k is the hopping parameter defined as k = 2(am +4) '^ith the lattice spacing a and quark mass iriq. We generate 
gauge ensembles at a fixed gauge coupling /? = = 6.2 with three different lattice sizes, xT = 24"^ x 48, 32'^ x 48 
and 48^ x 48. According to the Sommer scale [23, H^, P = 6.2 in the quenched approximation corresponds to a 
lattice cutoff of ~ 2.9 GeV, which may be marginal to handle low-lying cc mesons on the lattice. We compute 
the quark propagators at three values of the hopping parameter k = {0.1520,0.1506,0.1489}, which correspond to 
M^/Mp = 0.68,0.83,0.90. Kc = 0.1360 is reserved for the charm-quark mass. Simulation parameters and the number 
of the gauge samples are summarized in Tabled Preliminary accounts for L =24 and 32 cases are given in Ref. [2^ . 

For the update algorithm, we adopt the Metropolis algorithm with 20 hits at each link update. The first 10000 
sweeps are discarded for thermalization. The gauge ensembles in each simulation are separated by 1200 {L — 48), 800 
(L = 32) and 600 {L = 24) sweeps. For the matrix inversion, we use BiCGStab algorithm and adopt the convergence 
condition |7'| < lO^'* for the residues. We calculate smear-point quark propagators with a box-type smeared source 
which are located at tgrc — 6 for the light quarks and at tgrc = 5 for the charm quark. To enhance the coupling to 
ground states of hadrons, we choose a 24^ box source, of which spatial lattice size is about 1.6 fm, on all three volumes, 
instead of a wall {L^ box). Since the box- type smeared source is gauge variant. Coulomb gauge fixing is carried out by 
a combination of an SU(2) subgroup method and an overrelaxed steepest descent method. Details of this procedure 
may be found in Ref. [30|. To perform the precise parity projection for the nucleon, we adopt a procedure to take 
an average of two quark p ropagators which are subject to periodic and antiperiodic boundary conditions in time as 
described in Refs. 



A. Energy shift AE 

Hadron masses computed with the conventional single exponential fit are summarized in Tabic Hll Kc = 0.1360 
reproduces the mass of J/tp{3Q97) approximately. The mass of pseudoscalar meson becomes Af^ w 0.6 GeV for 
K = 0.1520, 0.9 GeV for k = 0.1506 and 1.2 GeV for k = 0.1489 in our simulations. 

The energy shifts of two-hadron systems {J/ip-n, J/ip-p and J/tJj-N) are determined by the large-i behavior of the 
ratio Rj/t(, h defined in Eq. I|10|l . To find appropriate temporal windows for fitting, we define the effective energy shift 

by 

Ai^effW-ln/^^^^, (18) 

^J/ijj-hV' + ij 

which should show a plateau for large Euclidean time (t ^ tgi-c) and approaches to the true energy shift AE, if the 
statistics is sufficiently large. 

As typical examples of our simulations, in Fig. ^ we show the effective mass AEcB{t) as a function of t for L — 24 
(top panels), 32 (middle panels) and 48 (bottom panels) at k = 0.1506. In the figure, the left panels are for the J/tp-n 
channel, the middle ones for the spin-0 J/ip-p channel, and the right ones for the spin- 1/2 J/tp-N channel. Statistical 
uncertainty of AEcs{t) is estimated by a single elimination jack-knife method. 

As seen from this figure, the effective energy shift is quite small compared to the mass of each hadron. Nevertheless, 
strong correlation between the numerator and denominator in the ratio Rj/^.h can help to expose such a tiny shift. 
An important observation is that the negative energy shift is found in all channels of Fig.^ which implies that there 
exist attractive interactions between J ftp and light hadrons (tt, p and N). 

One should, however, notice that the signals become quite noisy for t far away from the sources at tgrc — 5 and 
6. Therefore, we need to make appropriate choice of the fitting window in [tmin, imax] to extract AE from the data. 
Since the two-point correlation functions are well dominated by ground state hadrons { J/ip, tt, p and N) for t > 20, 
a reasonable choice of the lower bound is imin ~ 20. The upper bound may be chosen so that we have a reasonable 
value of x^/dof in the fitting interval. Three horizontal solid lines in each panel of Fig. ^represent the fitting range 
thus determined together with the value of AE and its la deviation obtained by a covariant single exponential fit to 
the ratio Rj/^.u- The fitting range is chosen to be the same for the data with equal L. The results of energy shift 
obtained in the above procedure are summarized in Tables UlTl II VI and IVl 



B. Quark mass dependence of the energy shifts 



Fig.|3shows the quark mass dependences of the energy shift AE for L = 24 (upper panels), L = 32 (middle panels) 
and L = 48 (lower panels). The left panels are for the J/tp-ir channel, the middle ones for the spin-0 J/t/j-p channel 
and the left ones for the spin-1/2 J/ip-N channel. Open circles show the results evaluated at three values of the 
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hopping parameter for the hght hadrons, k = 0.1520, 0.1506 and 0.1489. We do not find appreciable spin dependence 
of the energy shift within the error bars in the J/ip-p and J/ip-N channels. 

The J/iJj-TT channel has a special feature in the sense that the low-energy pion decouples from any hadrons in the 
chiral limit as discussed in Sec. 1. In other words, AE of the J/tp-i: system should vanish in simultaneous chiral and 
infinite volume limits. The top, middle and bottom panels for the J/ip-ir channel in Fig. |2] indeed suggest such a 
tendency. This will be quantified later in Sec. IIII Dl 

To make an extrapolation of AE to the chiral limit, we adopt a simple formula: 

aAE = c + c'iaAU)"^, (19) 

where c and c' are determined numerically from the data. In Fig. [2 full circles represent the energy shifts linearly 
extrapolated by Eq. (|19|) with the physical pion mass squared (M^r =140 MeV). For the J/ip-p and J/ip-N channels, 
we show also the results (full squares) obtained from the weighted average of the data for two heavy-quark masses 
(k = 0.1506 and 0.1489). The purpose is to check the sensitivity from the lightest data points which have relatively 
large error bars. We find that results of full circles and squares are in agreement with each other within the errors, 
and hence we quote the results evaluated by Eq. l{ll?)l in the following sections. 

C. Volume dependence of the energy shifts 

In Fig. 121 we show the volume dependence of the energy shift AE' at the physical point (M7r=140 MeV). The left 
panel is for the J/t/i-tt channel, the middle for the spin-0 J/i^-p channel and the right for the spin-1/2 J/4>-N channel. 
The horizontal axis denotes the spatial size L and the vertical axis the energy shift AE. The dashed lines indicate the 
lower boundary given in Eq. 1)15(1 : If the data points come below the dashed lines, the 1/L-expansion of the phase-shift 
formula is no longer justified. 

The volume dependence of the energy shift AE is a key quantity to discuss the validity of the large-L expansion of 
the phase-shift formula. We find that the absolute value of the energy shift in the J/ip-ir channel of Fig. O decreases 
monotonically from L = 24 to 48, which is consistent with the leading behavior shown in Eq. (|14|l . On the other 
hand, the energy shift in the J/ip-p and J/ip-N channels do not show such monotonic decrease due to rather small 
values of AE at L = 24. This implies that the lattice size L = 24 {La ~ 1.6 fm) is too small to satisfy the condition 
L > 2R ^ 2ao so that the convergence of the power series in Eq. (|14() is questionable or even the phase-shift formula 
in Eq. Ijl^l) itself is invalid. Thus, in the J/ip-p and J/ip-N channels, we do not take the results of L = 24 as physical 
in the following analyses. 

D. Volume dependence of scattering lengths 

In the left panel of Fig. 0] we show the volume dependence of the s-wave scattering length a^^^''^ at the physical 
point (M^ =140 MeV) obtained by applying the phase-shift formula H13() to AE in Fig.O The horizontal axis is the 
spatial size L and the vertical axis is the scattering length. By taking into account the error bars for L = 48 due to 
the limited statistics, we do not find appreciable L dependence for the J/i/^-tt scattering lengths. Also, the scattering 
length is positive, which implies attraction between J/iJj and tt at low energy. 

In the middle panel of Fig. ^ we show the volume dependence of the s-wave scattering length Qq^^''' obtained by 
Eq. (|13() . Within the error bars, different spin states have similar scattering lengths. The statistics in the L = 24 
case is much better than L = 32 and 48. However, the L = 24 data is contaminated by significant finite size effect as 
discussed in Sec. IIII d on the basis of Fig.O The scattering lengths for L — 32 and 48 show the positive values. 

In the right panel of Fig. 0] we show the volume dependence of the J/ip-N scattering length Qq^^'^ obtained by 
Eq. H13|l . Similar to the case in the J/ip-p channel, we do not find appreciable spin dependence within the error bars. 
Here the same discussion as the J/ip-p channel applies to the L = 24 data. Although we need to increase statistics 
to make definite conclusion, the scattering length in J/ip-N channel for L = 32 and 48 is positive and is considerably 
larger than that in the J/ifj-iT channel. 

In Table IVTI we summarize the s-wave scattering lengths in lattice units at the physical point (Af7r=140 MeV) and 
at the chiral limit (M^ = 0). The table shows a clear tendency, 

< a^/'^-" « 4""-" < ai'^-"". (20) 

This indicates that (i) there is always attraction between J /-ip and light hadrons, and (ii) there may be some relation 
between the strength of the attraction and the light-hadron properties such as the spatial size of light hadrons, or the 
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number of constituent quarks in the light hadrons. To clarify the point (ii) further, it will be necessary not only to 
make high statistics simulations but also to carry out systematic calculations for the interactions of J with various 
light hadrons. It is also worth mentioning here that the s-wave J/^'-tt scattering length obtained at each volume 
is consistent with zero in the chiral limit within 1 or 2 standard deviations. This is in accordance with the current 
algebra result discussed in Eq. 

E. Scattering lengths in the physical units 

In Table Ivnl the results of the s-wave scattering lengths, aj^~^ ^ for each h with different spin states are tabulated, 
where "SAV" implies the results after the spin average. We show the results of two different methods to extract the 
scattering length from the energy shift. 

"PSF (=Phase Shift Formula)" implies the method by using the phase-shift formula in Eq. H13|) without making 
the 1/L-expansion. This is what we utilized to extract the numbers in Table IVll In this case, the final results of 
the scattering lengths are obtained by an average over the scattering lengths obtained in different volumes. In an 
alternative method indicated by "LLE (^Leading large-L Expansion)" we take the leading 1/L formula in Eq. I|14|l . 
AE{L) = A ■ L~^, and calculate A by fitting the volume dependence of the energy shift. The scattering length is 

then extracted as a^^'^ = —[lAjl-K. The values of A obtained through this method are summarized in Table IVlTTl 
Theoretically, PSF is more rigorous than LLE. Nevertheless, the latter is useful to estimate a systematic error caused 
by a finite L. 

As we have discussed in Sections IIII CI and IIII Dl we know that L = 24 is not large enough to extract any useful 
information in the J/ip-p and J/tp-N channels. Therefore, in both PSF and LLE, we use only the data for L = 32 
and 48 in the J/ip-p and J/ip-N channels, while we use all data {L = 24,32 and 48) in the J/tp-ir channel. 

As seen in Table fVIll the results obtained from two analysis (PSF and LLE) agree with each other within the errors. 
Especially, in the J/V'-tt case, both estimates are fairly consistent with each other. The table also shows the elastic 
cross sections at the threshold given by CTci = Attuq. We find that the s-wave scattering length in the J/tp-ir channel 
is quite small, 0.0119(39) fm, but definitely positive. Therefore, we conclude that the J /ip-i^ interaction is attractive 
at least at low energy. The soft pion theorem, where implies that the pion should decouple from any other hadrons 
in the chiral limit {M-^ = 0), is an account for such smallness of the s-wave scattering length. The J/ip-i: elastic cross 

section at the threshold is evaluated as o'H^''" = 0.018to oio This is consistent with a phenomenological analysis 
based on short distant QCD by Fujii and Kharzeev jl^ . 

Although statistical errors in J/tp-p and J/ip-N channels are quite large, we find that s-wave scattering lengths in 
these channels are likely to be positive and are at least an order of magnitude larger than the pion case. As for the 
J/ip-N channeljthe scattering length from QCD sum rules by Havashigaki 12.1, gluonic van der Waals interaction by 
Brodsky et al. and the QCD multipole expansion by Voloshin et al. [lj| show that the s-wave scattering length 
is about 0.1 fm, 0.25 fm, and 0.37 fm or larger, respectively. Our central value of the scattering length 0.71 fm from 
PSF (0.39 fm from LLE) is comparable or even larger than those estimates, but we need to increase statistics of our 
data to draw solid comparison. 

F. Possible contaminations from channel mixings 

Before closing the discussion on J/i/i-hadron scattering lengths, we comment on possible contaminations from the 
r]c-h and DD states. If there were open channels with lower energy than the J/f/i-hadron system at the threshold, the 
Liischer's formula for extracting the scattering phase shift from AE in the desired channel is not applicable. 

The J/ip-TT system is free from the contamination of the tJc-tt subthreshold state. This is because the s-wave tJc-tt 
state has different total spin from the s-wave J/ip-n state. On the other hand, for the J/ip-p case, we cannot exclude 
the possible contamination of the ijc-p state to the spin-1 J/tp-p state and that of the DD state to the spin-0 J/ip-p 
state. Only the highest spin state (spin-2) is definitely free from such contaminations. This is true also for the case 
of J/ip-N system: The highest spin state (spin-3/2) is free from the contamination of the rjc-N state. 

Thus, strictly speaking, the spin-1 J/tp-ir, the spin-2 J/tp-p and the spin-3/2 J/tp-N states are the safe channels 
in determining the s-wave scattering phase shifts from the Liischer's formula, although in our simulations we do not 
find any appreciable difference among different spin channels within the error bars. 
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G. r^c-hadron interactions 

The scattering length of rjc with Ught hadrons can be performed in exactly the same way as the J/^ case by using 
the conventional rjc interpolating operator as in Sec. m 

O'^^ix^t) =Ca{x,t)-f5Ca{x,t). (21) 

In the present paper, we neglect the disconnected diagrams in both the two-point function of rj^ and the related four- 
point functions as we have done in the case of J/'0- From the theoretical point of view, neglecting the disconnected 
diagrams in the tJc channel is less justified than that in the J/tp channel. On the other hand, it is numerically 
observed that the contributions from the disconnected diagrams are negligibly small for the charm quark even in the 
pseudoscalar channel [25ll2^ . 

In our simulations the mass difference between the J/tp and 7]c is about 30 MeV, which is considerably smaller than 
the experimental value, 116 MeV. It is known that the hyperfine split ting is sensitive to the quenched approximation 
and the leading discretization error of the Wilson fermion action |25Ll2a |. Therefore, the spin-dependent observables 
in our simulations should be taken with caution. 

Under these reservations, the volume dependences of the energy shift and the scattering length are shown in Figs. [51 
and El respectively. Qualitative features of the figures are similar to those in the J/ ip case: The energy shift is 
definitely negative in all three channels at L = 32 and 48. Also, we find < clq" '^ <C a^" '' < Uq"' . In the rjc-p and 
r]c-N channels, considerable finite L effect may exist at L — 24. All analyses are made in the same way as the case of 
J/i(j. Table Hxl contains two types of the chirally extrapolated value of the s-wave scattering length in all channels for 
each volume. In Table 1x1 we compare results obtained from two different analyses (PSF and LLE). The coefficients 
A in LLE are tabulated in Table Ixll 

The s-wave tJc-tt scattering length turns out to be 0.0113±0.0035 fm and the corresponding elastic cross section at 
the threshold is cr^f^ = 0.016tooo8 both of which are very close to the values in the J/ip-ir channel. Although 
we have large error bars, the scattering lengths in the rjc-p and rjc-N channels are also comparable with those in the 
J/ip-p and J/ip-N channels. 

IV. SUMMARY AND CONCLUDING REMARKS 

In this paper, we have studied the interactions of J/ip and rjc with light-hadrons at low energy in quenched lattice 
QCD simulations. We calculated the scattering lengths of J/ip-ir, p, N and ric-7T,p,N from the energy shifts AE of 
two hadrons in a finite periodic box L'^ by using the Liischer's phase-shift formula. We employed the lattice spacing 
a ~ (2.9 GeV)~^, three light-quark masses M^^/Mp ~ 0.68, 0.83, 0.90, and three lattice sizes La « 1.6, 2.2 and 3.2 fm. 
Their values were utilized to perform analyses on quark-mass and lattice-size dependences of the energy shifts and 
the scattering lengths. 

We did not find appreciable quark-mass dependence of AE in all channels, while in the J/ip-ir channel it turned 
out that there is a tendency that the energy shift approaches zero in simultaneous chiral and large-volume limits as 
indicated in the current algebra analysis. 

As for the volume dependence of the energy shifts, the behavior, which is expected from the asymptotic 

expansion of the Liischer's formula in terms of l/L, was seen in the J/ip-n and rjc-TT channels. On the other hand, in 
the J/ifj-p, N and i]c-p, N channels, the data at L = 24 do not follow the behavior [33|. Indeed, La = 1.6 fm is 
not large enough to accommodate even a single hadron such as N 34] . 

On the basis of the above observation, we adopted the data for L = 24, 32 and 48 in the J/ip{r]c)-Tr channel and the 
data for L = 32 and 48 in the J/tp{ric)-p and J/4>{ric)-N channels. Then we applied the Liischer's formula without 
1/L-expansion (which we call PSF) and with 1/L-expansion (which we call LLE). The resultant s-wave scattering 
lengths obtained from AE show that (i) the interaction of charmonia with light hadrons is always attractive, (ii) 
the scattering lengths in the J/il){r]c)-TT channels are quite small compared to other channels as consistent with the 
soft-pion theorem, and (iii) the scattering lengths in the J/'ip(ric)-p and J/ip{r]c)-N channels are at least an order of 
magnitude larger than those in the J/'il){ric)-TT channels. 

The s-wave J/ip-ir {rjc-Tr) scattering length is determined as 0.0119±0.0039 fm (0.0113±0.0035 fm) and the cor- 
responding elastic cross section at the threshold becomes 0.018^^010 ™b (0.016lo;oJ8 mb). On the other hand, the 
J/ip-N {rjc-N) spin-averaged scattering length is 0.71±0.48 fm (0.70±0.66 fm) which has still a large statistical errors. 
Nevertheless, our result in the J/tp-N channel may be compared with estimates from QCD sum rules (~ 0.1 fm), 
from the gluonic van der Waals interaction 0.25 fm) and from the QCD multipole expansion (~ 0.37 fm or larger). 

To have more quantitative understanding of the scattering lengths, in particular, of their magnitudes and spin 
dependences, we need to accumulate more statistics. Since our data for the J/ip{r]c)-p{N) channels indicate that 
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the scattering length is rather large so that the large-L expansion is not fully justified even for L = 32 and 48. 
Therefore, it is important to go to larger volume to reduce the systematic errors due to the finite volume. Simulations 
with dynamical quarks are also an important future direction to be explored for full knowledge of the charmonium 
interactions. 
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Appendix A: Spin projection 

In the cases of the s-wave J /ip-p and J/tp-N scatterings, there are different spin states: Spin-0, 1 and 2 states for 
the J/ip-p and spin-1/2 and 3/2 states for the J/ip-N . Therefore, the appropriate spin projections are required to 
disentangle each spin contribution from four-point correlation functions. 

For the J/i/'-p system, the four point function, Eq. can be expressed by the orthogonal sum of spin-0, spin-1 
and spin-2 components: 

Gilti'^ = G'it)P?rM + G\t)Pl^.M + G'{t)P^,ki (22) 
with spin projection operators defined by 

(23) 
(24) 
(25) 



pO 


— -^SijSkl, 






= ^{StkSji 




p2 


= ^{StkSji 


+ ^il^jk) ~ -^SijSkl 



Respective spin-projected correlators are given as 

3 

3 



- {GiS'i^-GiS'it}), (27) 

- ToT (GilS^'it) + GiS^'it) - iGiltr^iti^ , (28) 



where indices i and j should be summed over all spatial directions. 

The four point function for the J/tp-N system, Eq. (O, can be also decomposed into spin-1/2 and spin-3/2 compo- 
nents as 



Here, spin projection operators for spin-1/2 and spin-3/2 are given by 



Gi'^-''{t) = G^I\t)P]l^ + G^I\t)P^I\ (29) 



= ^7.7,, (30) 
pT = %-^7.7,- (31) 
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Then, each spin part can be projected out as 

G'/'it) - E i^y^G/^"(t) = I E 7.7.G;/'^-"(0, (32) 

G^t) = \ E P^/'G-H^'-^'W - I E - I E 7.T,G^"-"(0, (33) 

i.j — l i—1 — ^ 

where indices i and j are also summed over ah spatial directions. Recall that respective contributions G^^^{t) and 
G'^^^{t) possess non-trivial Dirac structure. To extract the particle contribution, we need to take a trace with the 
projection operator (1 +74)72 |3ll |. 
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TABLE I: Simulation parameters in this study. The Sommer parameter ro = 0.5 fm is used to fix the scale '57, "2^. 





a [fm] 


[GeV] 


Lattice size {I? x T) 


~ La [fm] 


Statistics 


6.2 


0.06775 


2.913 


24'' X 48 


1.6 


161 








32^ X 48 


2.2 


169 








48^ X 48 


3.2 


53 



TABLE IL Fitted masses of pseudoscalar, vector and nucleon states in lattice units. We perform a covariant single exponential 
fit to two-point functions of each hadron in respective fitting ranges. The vector-meson mass at k =0.1360 in the physical unit 
is close to the mass of J/'(/'(3097) and hence the value is reserved for a charm quark. Other three hopping parameters are for 
light quarks. The box-to-point quark propagators are used in the present study, while the point-to-point quark propagators 
are used in Ref. '55] with different gauge configurations. The hadron masses in two approaches agree well with each other. 



X T 


Fitting range 


ft 




aMp 


aMN 


24^ X 48 


[22,31] 


0.1360 
0.1489 
0.1506 
0.1520 


1.017(1) 
0.415(1) 
0.314(2) 
0.213(2) 


1.027(1) 
0.461(2) 
0.383(3) 
0.322(6) 


1.584(2) 
0.722(4) 
0.589(5) 
0.467(10) 


32^ X 48 


[25, 34] 


0.1360 
0.1489 
0.1506 
0.1520 


1.020(1) 
0.416(1) 
0.315(1) 
0.213(1) 


1.030(1) 
0.463(1) 
0.383(2) 
0.317(4) 


1.594(2) 
0.723(4) 
0.592(4) 
0.473(8) 


48^ X 48 


[17,26] 


0.1360 
0.1489 
0.1506 
0.1520 


1.016(1) 
0.415(1) 
0.314(1) 
0.212(1) 


1.026(1) 
0.460(2) 
0.380(2) 
0.311(3) 


1.581(3) 
0.719(4) 
0.586(4) 
0.465(5) 


TABLE IIL Energy shifts Ai5 in all J/ip-h channels on the lattice with L — 24. The energy shifts are obtained from the ratios 
Rj/ii-h by the single exponential fit. The fitting range is chosen to be 22 < t < 31 for all channels so as to be the same as that 
employed for spectroscopy of single hadrons for L = 24. The labels "phys." and "chiral" show the values extrapolated to the 
physical point {AItt = 140 MeV) and to the chiral limit, respectively, with the linear quark-mass dependence in Eq. 119|l. 


ft 












0.1489 
0.1506 
0.1520 


-0.0013(4) 
-0.0014(4) 
-0.0016(5) 


-0.0020(4) -0.0016(4) 
-0.0016(6) -0.0012(6) 
-0.0006(13) 0.0003(13) 


-0.0012(4) 
-0.0007(6) 
0.00010(13) 


-0.0028(9) 
-0.0023(10) 
-0.0016(16) 


-0.0032(9) 
-0.0028(11) 
-0.0023(17) 


phys. 
chiral 


-0.0016(5) 
-0.0016(5) 


-0.0008(12) 0.0000(12) 
-0.0008(12) 0.0001(12) 


0.0007(12) 
0.0007(12) 


-0.0014(17) 
-0.0014(17) 


-0.0022(17) 
-0.0022(18) 
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TABLE IV: Energy shifts for all J/tp-h channels on the lattice with L — 32. The fitting range is chosen to be 25 < t < 34 for 
all channels. The labels are the same as in Table ITTTl 

















0.1489 
0.1506 
0.1520 


-0.0010(3) 
-0.0009(3) 
-0.0008(3) 


-0.0021(5) 
-0.0023(6) 
-0.0027(11) 


-0.0019(5) 
-0.0021(6) 
-0.0023(11) 


-0.0016(5) 
-0.0017(6) 
-0.0018(10) 


-0.0032(15) 
-0.0028(14) 
-0.0029(19) 


-0.0035(16) 
-0.0033(16) 
-0.0039(22) 


phys. 
chiral 


-0.0007(3) 
-0.0007(3) 


-0.0028(11) 
-0.0028(11) 


-0.0024(10) 
-0.0024(10) 


-0.0019(10) 
-0.0020(10) 


-0.0027(19) 
-0.0027(19) 


-0.0036(21) 
-0.0036(22) 


TABLE V: Energy shifts for all J/ip-h channels on the lattice with L 
all channels. The labels are the same as in Table lUll 


— 48. The fitting range is chosen to be 17 < t < 26 for 


K 


aAEj/^_^ 












0.1489 
0.1506 
0.1520 


-0.0010(3) 
-0.0008(2) 
-0.0005(2) 


-0.0015(4) 
-0.0014(3) 
-0.0011(3) 


-0.0014(4) 
-0.0013(3) 
-0.0009(3) 


-0.0013(4) 
-0.0012(3) 
-0.0009(3) 


-0.0025(8) 
-0.0020(6) 
-0.0017(6) 


-0.0026(8) 
-0.0022(6) 
-0.0018(6) 


phys. 
chiral 


-0.0004(2) 
-0.0004(2) 


-0.0010(4) 
-0.0010(4) 


-0.0008(4) 
-0.0008(4) 


-0.0008(4) 
-0.0008(4) 


-0.0014(7) 
-0.0014(7) 


-0.0015(7) 
-0.0015(7) 



TABLE VL Results of the chirally extrapolated values of s-wave scattering lengths in lattice units for all J/tp-h channels. The 
labels are the same as in Table HTn 







j/ii-jT 


J/i'-p 

"0 










xT 






spin-0 


spin-1 


spin- 2 


spin- 1/2 


spin-3/2 


24^ X 48 


phys. 
chiral 


0.16(5) 
-0.0029(18) 


0.38(62) 
0.38(62) 


-0.03(56) 
-0.03(56) 


-0.34(53) 
-0.34(53) 


0.97(1.32) 
0.96(1.31) 


1.66(1.59) 
1.65(1.57) 


32'' X 48 


phys. 
chiral 


0.17(7) 
-0.0013(10) 


4.4(2.4) 
4.4(2.4) 


3.6(2.0) 
3.6(2.0) 


2.8(1.7) 
2.7(1.6) 


6.6(7.5) 
6.5(7.4) 


11.7(13.4) 
11.5(13.0) 


48^ X 48 


phys. 
chiral 


0.33(17) 
-0.0014(15) 


5.1(2.4) 
5.0(2.4) 


3.8(2.1) 
3.7(2.0) 


3.6(2.1) 
3.6(2.1) 


12.5(11.2) 
12.3(11.0) 


14.0(12.8) 
13.8(12.5) 



TABLE VIL The s-wave scattering lengths of JZ-i/i-hadrons and the elastic cross sections at the threshold in the physical unit. 
For the J/ip-n channel, the data for all three volumes of L = 24, 32 and 48 are used, while for the J/ip-p a-nd J/tp-N channels, 
only the data for L = 32 and 48 are used. In the table, "PSF" and "LLE" stand for the Phase Shift Formula and Leading 
large-L Expansion, respectively. "SAV" stands for the spin-averaged value, |[5(ao)2 + 3(ao)i -I- (ao)o], for the J/ip-p channel 
and i[2(ao)3/2 -I- (10)1/2] for the J/ip-N channel. 







From PSF 




From LLE 




Channel 


Spin 


fflo [fm] 


(Jci [mb] 


do [fm] 


(Joi [mb] 


J/tp-n 


1 


0.0119±0.0039 


0.01811!;;;;^, 


0.0119±0.0025 


U.Ui»_o.oo7 


J/i'-p 





0.32±0.12 


12.9+^^e" 


0.23±0.06 


6-6_3.2 




1 


0.25±0.10 




0.19±0.06 


4 6+^-^ 

^•"-2.4 




2 


0.21±0.09 




0.17±0.06 


r-l-2.8 
'^■'^-2.0 




SAV 


0.23±0.08 




0.18±0.05 


4 1+2-9 


J/iP-N 


1/2 


0.57±0.42 


41 + 38 


0.35±0.15 






3/2 


0.88±0.63 




0.43±0.16 






SAV 


0.71±0.48 




0.39±0.14 





TABLE VIII: Fitting parameter A as the coefficient of the leading term in the large-L expansion of the energy shift. For the 
J/tp-n channel, we used all energy shifts at L = 24, 32 and 48 to extract the fitting parameter A. On the other hand, for the 
J/i'-p and J/ip-N channels, we excluded the energy shifts at L = 24 for evaluating the parameters A. 









41 

^J/i'-p 


42 


4I/2 
^.//V'-iV 


43/2 


phys. 


-24(5) 


-101(28) 


-84(26) 


-73(25) 


-115(48) 


-140(52) 
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TABLE IX: Results of the chiral extrapolated values of s-wave scattering lengths in lattice units for all jyc-hadron channels. 
The labels are the same as in Table ITTTl 



X T 


"o 


"o 


Vc-N 


24-' X 48 


phys. 
chiral 


0.15(5) 
-0.0025(15) 


-0.25(48) 
-0.25(48) 


1.10(1.29) 
1.09(1.28) 


32'' X 48 


phys. 
chiral 


0.15(5) 
-0.0013(9) 


2.8(1.6) 
2.7(1.6) 


8.8(9.2) 
8.7(9.0) 


48'' X 48 


phys. 
chiral 


0.33(17) 
-0.0014(14) 


3.4(1.8) 
3.4(1.8) 


12.4(10.5) 
12.2(10.2) 



TABLE X: The s-wave scattering lengths of r/c-hadrons and the elastic cross sections at the threshold in the physical unit. For 
the rjc-n channel, the data for all three volumes of L = 24, 32 and 48 are used, while for the rjc-p and ric-N channels, only the 
data for i = 32 and 48 are used. The labels are the same as in Table fvTTl 







From PSF 




From LLE 




Channel 


Spin 


ao [fm] 


(Jci [mb] 


ao [fm] 


(Jci [mb] 







0.0113±0.0035 


016+" "^^ 


0.0112±0.0024 


016+"-""** 


Vc-P 


1 


0.21±0.11 




0.16±0.05 


■^•^-1.8 


Vc-N 


1/2 


0.70±0.66 


62±^^^ 


0.39±0.14 





TABLE XL Fitting parameter A as the coefficient of the leading term in the large-L expansion of the energy shift. For the 
rjc-TT channel, we use all energy shifts at L = 24, 32 and 48 to extract the fitting parameter A. On the other hand, for the rjc-p 
and rjc-N channels, we exclude the energy shifts at L = 24 for evaluating the parameters A. 







Aric-p 


Arj^-N 


phys. 


-23(5) 


-72(23) 


-127(47) 
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0.000 




-0.008 



10 20 30 
t : time slice 



0.000 



-0.004 



-0.008 




0.000 



-0.008 




-0.002 



-0.006 




10 20 30 
t : time slice 



40 50 




L = 32, Spin 
o K = 0.1506 



10 20 30 40 
t : time slice 



-0.006 




20 30 40 50 
t : time slice 



-0.002 



-0.006 




10 20 30 
t : time slice 



-0.002 



-0.006 




L = 48, Spin 

o K = 0.1 506 



10 20 



40 50 



-0.006 




L = 48, Spin 1/2 
o K = 0.1 506 



10 20 



40 50 



t : time slice 



t : time slice 



FIG. 1: The effective energy shifts A.E in lattice units as a function of the time slice t. Hopping parameters of light hadrons 
and the J/tp are fixed to be k = 0.1506 and 0.1360, respectively. Top (middle, bottom) panels are for L =24 (32, 48). Left 
(middle, right) panels are for the J/^J-tt {J/ilj-p in spin-0, J/ip-N in spin-1/2) channels. Source locations of light hadrons and 
the J/'ip are at tsrc = 6 and 5, respectively. 



15 



0.002 rr 



-0.008 




0.05 0.10 0.15 
: Pion Mass Squared 



0.002 n- 



-0.006 



-0.008 



L = 24, Spin 
o Data 
• Linear 
■ Ave. 
I 



0.05 0.10 0.15 
m/ : Pion Mass Squared 




-0.008 



0.05 0.10 0.15 0.20 
m/ : Pion Mass Squared 



-0.004 



-0.006 



1 1 1 1 1 








- L = 32 






o Data 






• Linear 
1 1 


1 





0.00 0.05 0.10 0.15 
: Pion Mass Squared 



5) -0.002 



-0.006 



L = 32, Spin 
o Data 
• Linear 
■ Ave. 
I 



0.00 



0.05 0.10 0.15 
: Pion Mass Squared 



B> -0.002 



-0.006 



L = 32, Spin 1/2 
o Data 
• Linear 
■ Ave. 
i_ 



0.05 0.10 0.15 0.20 
: Pion Mass Squared 



0.000 



-0.002 



-0.004 



1 1 1 1 1 






- L = 48 




o Data 




• Linear 
1 1 


1 1 1 



0.05 0.10 0.15 
m^^ : Pion Mass Squared 



-0.004 



1 1 1 1 1 






L = 48, Spin 




- o Data 




• Linear 




■ Ave. 

1 1 1 


1 1 



0.05 0.10 0.15 
niji^ : Pion Mass Squared 



-0.004 




0.05 0.10 0.15 0.20 
: Pion Mass Squared 



FIG. 2: The energy shifts AE as a function of the squared pion mass in lattice units. Upper (middle, lower) panels are 
for L = 24 (32, 48). Left (middle, right) panels are for the J/t/j-tt (J/f/j-p in spin-0, J/^-N in spin-1/2) channels. Open circles 
are the values fitted by the criterion obtained from the four-point function at the hopping parameters of light hadrons (k — 
0.1489, 0.1506 and 0.1520). FuU circles are the values extrapolated to the physical point {AU =140 MeV) by using Eq. itT^ . 
Full squares are the values from weighted averages of heavy-quark results (k — 0.1489 and 0.1506). 
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FIG. 3: The energy shifts AE as a function of spatial size L in lattice units. Left (middle, right) panel is for the J/^-n (J/tp-p 
in spin-0, J/tp-N in spin-1/2) channel. Full and open circles represent the values obtained from the linear chiral-extrapolation 
to the physical point, and the weighted average, respectively. The dashed curves show the lower boundary for the convergence 
of the large-!/ expansion of AE. 
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FIG. 4: The scattering lengths as a function of the spatial size L in lattice units for physical pion mass (M,r = 140 MeV). Left 
(middle, right) panel is for the J/tp-n {J/ip-p, J/ip-N) channel. 
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FIG. 5: The energy shifts AE as a function of spatial size L in lattice units. Left (middle, right) panel is for the t/c-tt (ijc-p, 
rjc-N) channel. Full and open circles represent the values obtained from the linear chiral-extrapolation to the physical point, and 
the weighted average, respectively. The dashed curves show the lower boundary for the convergence of the large-L expansion 
of AE. 
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FIG. 6: The scattering lengths as a function of the spatial size L in lattice units for the physical pion mass {Mt, = 140 MeV). 
Left (middle, right) panel is for the r/c-vr {rjc-p, rjc-N) channel. 



